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Abstract 

In this paper we propose a collocation method for solving some well-known classes of Lane- 
Emden type equations which are nonlinear ordinary differential equations on the semi-infinite 
domain. They are categorized as singular initial value problems. The proposed approach is 
based on a Hermite function collocation (HFC) method. To illustrate the reliability of the 
method, some special cases of the equations are solved as test examples. The new method 
reduces the solution of a problem to the solution of a system of algebraic equations. Hermite 
functions have prefect properties that make them useful to achieve this goal. We compare 
the present work with some well-known results and show that the new method is efficient 
and applicable. 
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1. Introduction 

Many problems in science and engineering arise in unbounded domains. Different spec- 
tral methods have been proposed for solving problems on unbounded domains. The most 
common method is through the use of polynomials that are orthogonal over unbounded 
domains, such as the Hermite spectral and the Laguerre spectral methods []]-|_|. 

Guo [i- II] proposed a method that proceeds by mapping the original problem in an 



unbounded domain to a problem in a bounded domain, and then using suitable Jacobi 
polynomials to approximate the resulting problems. 
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Another approach is replacing the infinite domain with [-L, L] and the semi-infinite 
interval with [0, L\ by choosing L, sufficiently large. This method is named as the domain 
truncation [12] . 

Another effective direct approach for solving such problems is based on rational approx- 
imations. Christov 13| and Boyd [l4j, |l5| developed some spectral methods on unbounded 
intervals by using mutually orthogonal systems of rational functions. Boyd 15| defined a 
new spectral basis, named rational Chebyshev functions on the semi-infinite interval, by 
mapping it to the Chebyshev polynomials. Guo et al. Q introduced a new set of rational 
Legendre functions which are mutually orthogonal in L 2 (0, +00). They applied a spectral 
scheme using the rational Legendre functions for solving the Korteweg-de Vries equation on 



the half line. Boyd et al. [17] applied pseudospectral methods on a semi-infinite interval 



and compared the rational Chebyshev, Laguerre and the mapped Fourier sine methods. 

Authors of 18-23[ applied the spectral method to solve the nonlinear ordinary differen- 
tial equations on semi-infinite intervals. Their approach is based on the rational Tau and 
collocation methods. 

Lane-Emden type equations are nonlinear ordinary differential equations on semi-infinite 
domain. They are categorized as singular initial value problems. These equations de- 
scribe the temperature variation of a spherical gas cloud under the mutual attraction of 
its molecules and subject to the laws of classical thermodynamics. The polytropic theory of 
stars essentially follows out of thermodynamic considerations, that deals with the issue of 
energy transport, through the transfer of material between different levels of the star. These 
equations are one of the basic equations in the theory of stellar structure and has been the 

focus of many studies jUH, SB 0, H, H3, E3 • 



We simply begin with the Poisson equation and the condition for hydrostatic equilibrium 



46, 55-57 
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where G is the gravitational constant, P is the pressure, M(r) is the mass of a star at a 



certain radius r, and p is the density, at a distance r from the center of a spherical star [57 



The combination of these equations yields the following equation, which as should be noted, 



is an equivalent form of the Poisson equation [46j, I55H57] : 



1 d /r 2 dP\ 



r 2 dr ^ p dr ) ^ 

From these equations one can obtain the Lane-Emden equation through the simple assump- 
tion that the pressure is simply related to the density, while remaining independent of the 
temperature. We already know that in the case of a degenerate electron gas, the pressure 



2 



and density are p ~ P$ , assuming that such a relation exists for other states of the star, we 
are led to consider a relation of the following form 571 ] : 



P = Kp 1+ ™, 

where K and m are constants, at this point it is important to note that m is the polytropic 
index which is related to the ratio of specific heats of the gas comprising the star. Based 
upon these assumptions we can insert this relation into our first equation for the hydrostatic 
equilibrium condition and from this equation 46|, 1551457] we have 



g(m+l) j, _! 

A m 



r 2 dr V dr I 



where the additional alteration to the expression for density has been inserted with A rep- 
resenting the central density of the star and y that of a related dimensionless quantity that 
are both related to p through the following relation 46, 55-5?J 



p = \y m . 

Additionally, if place this result into the Poisson equation, we obtain a differential equation 
for the mass, with a dependance upon the polytropic index m. Though the differential 
equation is seemingly difficult to solve, this problem can be partially alleviated by the 
introduction of an additional dimensionless variable x, given by the following: 



ax, 



K(m + 1) % j_ x 

A m 

4ttG 

Inserting these relations into our previous equations we obtain the famous form of the Lane- 
Emden equations, given in the following: 



1 d 

x 2 dx 



x 



dx 



Taking these simple relations we will have the standard Lane-Emden equation with g(y) 
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y" + -y' + y r 

x 



0. 



x > 0. 



'1.11 



At this point it is also important to introduce the boundary conditions which are based 
upon the following boundary conditions for hydrostatic equilibrium and normalization con- 
sideration of the newly introduced quantities x and y. What follows for r = is 



0^x = 0, p = \^y(0) = l. 



;i.2) 
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As a result an additional condition must be introduced in order to maintain the condition 
of Eq. (11.21) simultaneously: 

2/(0) = 0. 

In other words, the boundary conditions follows 
y(0) = l, y'(0) = 0. 

The values of m which are physically interesting, lie in the interval [0,5]. The main difficulty 
in the analysis of this type of equation is the singularity behavior occurring at x = 0. 
Exact soloutions for Eq. f 1 1.1 1) are known only for m = 0, 1 and 5. For other values of m the 
standard Lane-Emden equation is to be integrated numerically. Thus we decided to present 
a new and efficient technique to solve it numerically. 

This paper is arranged as follows: 
In Section [2] we survey several methods that have been used to solve Lane-Emden type 
equations. In Section [31 the properties of Hermite functions and the way to construct the 
collocation technique for this type of equation are described. In Section 0] the proposed 
method is applied to some types of Lane-Emden equations, and a comparison is made with 
the existing analytic or exact solutions that were reported in other published works in the 
literature. Finally we give a brief conclusion in the last section. 



2. Methods have been proposed to solve Lane-Emden equations 



Recently, many analytical methods have been used to solve Lane-Emden equations, the 
main difficulty arises in the singularity of the equations at x = 0. Currently, most techniques 
which were used in handling the Lane-Emden-type problems are based on either series 
solutions or perturbation techniques. 

Bender et al. [24j proposed a new perturbation technique based on an artificial parameter 
5, the method is often called 5-method. 
Mandelzweig et al. 



251 ] used the quasilinearization approach to solve the standard Lane- 



Emden equation. This method approximates the solution of a nonlinear differential equation 
by treating the nonlinear terms as a perturbation about the linear ones, and unlike pertur- 
bation theories is not based on the existence of some small parameters. 

Shawagfeh 26| applied a nonperturbative approximate analytical solution for the Lane- 
Emden equation using the Adomian decomposition method. His solution was in the form of 
a power series. He used Pade approximants method 27], |28j| to accelerate the convergence 
of the power series. 

In [29( , Wazwaz employed the Adomian decomposition method 30|, 131 



framework designed to overcome the difficulty of the singular point. 
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with an alternate 
It was applied to 
used the modified 



the differential equations of Lane-Emden type. Further author of 
decomposition method for solving the analytical treatment of nonlinear differential equations 
such as the Lane-Emden equation 
Liao 
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provided an analytical algorithm for Lane-Emden type equations. This algo- 
rithm logically contains the well-known Adomian decomposition method. Different from all 
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other analytical techniques, this algorithm itself provides us with a convenient way to adjust 
convergence reg ions even without Pade technique. 

J.-H He 34( employed Ritz's method to obtain an analytical solution of the problem. By 
the semi-inverse method, a variational principle is obtained for the Lane-Emden equation. 



Parand et al. [19l . |2Q| . |23| presented two numerical techniques to solve higher ordinary 



differential equations such as Lane-Emden. Their approach was based on the rational Cheby 
shev and rational Legendre Tau methods. 



Ramos |35H38j solved Lane-Emden equations through different methods. Author of [36 



presented the linearization method for singular initial-value problems in second-order ordi- 
nary differential equations such as Lane-Emden. These methods result in linear constant- 
coefficients ordinary differential equations which can be integrated analytically, thus yielding 
piecewise analytical solutions and globally smooth solutions. Later this author [38] developed 
piecewise-adaptive decomposition methods for the solution of nonlinear ordinary differential 
equations. In [37| , series solutions of the Lane-Emden type equation have been obtained by 
writing this equation as a Volterra integral equation and assuming that the nonlinearities are 
sufficiently differentiable. These series solutions have been obtained by either working with 
the original differential equation or transforming it into an ordinary differential equation 
that does not contain the first-order derivatives. Series solutions to the Lane-Emden type 
equation have also been obtained by working directly on the original differential equation 
or transforming it into a simpler one. 

Yousefi 39j presented a numerical method for solving the Lane-Emden equations. He 
converted Lane-Emden equations to integral equations, using integral operator, and then he 
applied Legendre wavelet approximations. 

Bataineh et al. 40] presented an algorithm based on homotopy analysis method (HAM) jU 
to obtain the approximate analytical solutions of the singular IVPs of the Emden-Fowler 
type equation. 

In [42J, Chowdhury et al. presented an algorithm based on the homotopy-perturbation 
method (HPM) [43-45] to solve singular IVPs of time-independent equations. 

Aslanov [46j introduced a further development in the Adomian decomposition method 
to overcome the difficulty at the singular point of non-homogeneous, linear and non-linear 
Lane-Emden-like equations. 

Dehghan and Shakeri 55| applied an exponential transformation to the Lane-Emden type 
equations to overcome the difficulty of a singular point at x = and solved the resulting 



nonsingular problem by the variational iteration method [58|, [59 



Yildirim et al. 47| presented approximate-exact solutions of a class of Lane-Emden type 
singular IVPs problems, by the variational iteration method. 



Marzban et al. 48 



used a method based upon hybrid function approximations. They 
used the properties of hybrid of block-pulse functions and Lagrange interpolating polyno- 
mials together for solving the nonlinear second-order initial value problems and the Lane- 
Emden equation. 

Recently, Singh et al. |49j provided an efficient analytic algorithm for Lane-Emden type 
equations using modified homotopy analysis method, also they used some well-known Lane- 
Emden type equations as test examples. 



We refer the interested reader to |5CJ, l5jj for analysis of the Lane-Emden equation based 
on the Lie symmetry approach. 



3. Hermite functions collocation method 

Spectral methods have been successfully applied in the approximation of boundary value 
problems defined in unbounded domains. For problems whose solutions are sufficiently 
smooth, they exhibit exponential rates of convergence/spectral accuracy. We can apply dif- 
ferent spectral methods that are used to solve problems in the semi-infinite domains. One of 
these approaches is using Laguerre and Hermite polynomials/functions [1, 52 - 54 ] . Guo 



[5| suggested the Laguerre- Galerkin method for the Burgers' equation and Benjamin-Bona- 
Mahony (BBM) equation on a semi-infinite interval. In [7[ Shen proposed spectral methods 
using Laguerre functions and analyzed the elliptic equations on regular unbounded domains. 
Siyyam [8] applied two numerical methods for solving differential equations using the La- 
guerre Tau method. Maday et al. [6] proposed a Laguerre type spectral method for solving 
partial differential equations. Funaro and Kavian [2[ considered some algorithms by using 
the Hermite functions. Recently Guo j3] developed the spectral method by using Hermite 
polynomials. However it is not easy to perform the quadratures in unbounded domains, 
which are used in the Hermite spectral approximations. So the Hermite pseudospectral 
method is more preferable in actual calculations. Guo (54( developed the Hermite pseu- 



dospectral method for the Burgers' equation on the whole line. Guo et al. [52j considered 



spectral and pseudospectral approximations using Hermite functions for partial differential 



equations (PDEs) on the whole line to approximate the Dirac equation. Bao and Shen [53 
proposed a generalized-Laguerre-Hermite pseudospectral method for computing symmetric 
and central vortex states in Bose-Einstein condensates (BECs) in three dimensions with 
cylindrical symmetry. 



Collocation method [60[ has become increasingly popular for solving differential equa- 
tions. Also they are very useful in providing highly accurate solutions to differential equa- 
tions. In this paper, we employ the Hermite functions collocation (denoted by HFC) method 
to solve some well-known singular forms of the Lane-Emden type initial value problems di- 
rectly. 

3.1. Properties of Hermite functions 

In this section, we detail the properties of the Hermite functions that will be used to 
construct the HFC method. First we note that the Hermite polynomials are generally not 



suitable in practice due to their wild asymptotic behavior at infinities [61 
Hermite polynomials can be written in direct formula as follows: 



TT i \ T(n+1) ^ . /- — — nir, 

HJx) ~ —r-, —e x 11 cos W2n + lx 

V ; r (n 2 + 1) V 2 ' 



n 



n/2 e x 12 cos( v / 2^TTx 



2/n , i nix. 



2 
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Hence, we shall consider the so called Hermite functions. The normalized Hermite functions 
of degree n is defined by 

1 



HJx) 



-.e~ x /2 H n (x), n>0,xe 



\f2 n n\ 

Clearly, {H n } is an orthogonal system in L 2 



+00 



H n (x)H m {x)dx = y/nS mn , 



where 5 nm is the Kronecker delta function. In contrast to the Hermite polynomials, the 
Hermite functions are well behaved with the decay property: 



\HJx) 



0. 



as x 



00, 



and the asymptotic formula with large n is 



nir , 



H n {x) ~ n 4 cos(v / 2n + lx — '-^-). 
The three-term recurrence relation of Hermite polynomials implies 



H n+1 (x) =x^l ~^H n {x) - y/^Hn-iix), n > 1, 
H (x) = e- x ' 2 / 2 , Hi(x) = V2xe- X ' 2 / 2 . 
Using the recurrence relation of Hermite polynomials and the above formula lead to 

H' n {x) = V2^H n -!(x) - xH n (x) 



n H n ^ 1 (x) - J^-^H n+1 (x). 



and this implies 



yj n(n—l)-w 



H' n (x)H'Jx)dx = { 



y/(n+l)(n+2)7r 
2 ' 



0, 



m = n — 2, 
m = n, 

m = n + 2, 
otherwise. 



Let us define 

P N ■= { u : u = e~ x2/2 v,Vv G P N }, 

where P N is the set of all Hermite polynomials of degree at most N. 

We now introduce the Gauss quadrature associated with the Hermite functions approach. 

Let {xj}^ =0 be the Hermite-Gauss nodes and define the weights 



(N + l)H 2 N {x 3 ) 



0<j<N. 



Then we have 



f N 

/ p(x)dx = y ^p(x j )w j , Wp G P 2N +i- 

j=0 



For a more detailed discussion of these early developments see [62|, |63 



3.2. Approximations by Hermite functions 
Let us define A := {x\ — oo < x < 00} and 

Hn = span{H (x), Hi(x), H N (x)}. 

The L 2 (A)-orthogonal projection £/v : L 2 (A) - 
v G L 2 (A), 



"Hat is a mapping in a way that for any 



< £ N v -v,<f>>= 0, 
or equivalently, 



G U 



N ■ 



N 



£ N v(x) = ^viHj 



1=0 



To obtain the convergence rate of Hermite functions we define the space H r A {h) defined by 

H r A (h) = {v\v is measurable on A and ||u|| r ^ < 00}, 

and equipped it with the norm \\v || rj 4 = ||A r ?;||. For any r > 0, the space H A (A) and its 
norm are defined by space interpolation. By induction, for any non-negative integer r we 
can write 



A r v{x) = ^(x 2 + l) {r - k)/2 p k {x)d k x v{x) 



k=0 



where pk{x) are certain rational functions which are bounded uniformly on A. Thus we have 
_^ t 1/2 

\v\\ r ,A < C ( J2 II + l) (r - k)/2 Pk(x)d k x V 



K k=0 



Theorem 3.1. For any v G H A (A), r > 1 and < \i < r, the following can be obtained: 



£NV-V\L< C ATV3+(/,-l)/2 



\r,A' 



(3.1) 



Proof. A complete proof is given by Guo et al. [52|. Also same theorems have been 



proved by Shen et al. [61 
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3. 3. Hermite functions transform 

As mentioned before, Lane-Emden type equations are defined on the interval (0, +00); 
but we know properties of Hermite functions are derived in the infinite domain (—00, +00). 
Also we know approximations can be constructed for infinite, semi-infinite and finite inter- 
vals. One of the approaches to construct approximations on the interval (0, +00) which is 
used in the current paper, is to use a mapping, that is a change of variable of the form 

lo = cj)(z) = ln(sinh(fe)), (3.2) 

where k is a constant. 

The basis functions on (0, +00) are taken to be the transformed Hermite functions, 

H n (x) = H n (x) o <p(x) = H n ((f)(x)), 
where H n (x) o <fi(x) is defined by H n (x)((j)(x)). The inverse map of 00 = <f)(z) is 

z = ^(co) = I We" + Ve 2 " + 1). (3.3) 
k 



Thus we may define the inverse images of the spaced nodes {xj} x '. = _ OQ as 
T = {0 _1 (t) : -00 < t < +00} = (0, +00), 

and 



7 ln(e^ + y/e 2x s + l), J =0,1,2,... 



Xj = <p [Xj) k 

Let w(x) denotes a non-negative, integrable, real-valued function over the interval T. We 
define 

L^(r) = {v : T — > R I t>is measurable and || v\\ w < 00}, 

where 

1 



|| v\\ w = yj I v{x) | 2 w(x)dx 
is the norm induced by the inner product of the space L 2 W (r), 

/>oo 

<u,v> w = / u(x)v(x)w(x)dx. (3.4) 



Thus {i? n (x)}„ e N denotes a system which is mutually orthogonal under (j3.4j) . i.e., 

< H n (x),H m (x) > W ( X )= y/n5 nm , 

where w(x) = coth(a;) and 5 nm is the Kronecker delta function. This system is complete in 
L 2 W (r). For any function / G L^(T) the following expansion holds 

9 



+N 
jfe=— JV 



with 

, _ < f(x),H k (x) > w(x .) 
11^(^)11^) 

Now we can define an orthogonal projection based on the transformed Hermite functions as 
given below: Let 

U N = span{H ld (x),H 1 (x), ...,H n (x)}. 

The L 2 (r)-orthogonal projection £n '■ L 2 (T) — > "Hat is a mapping in a way that for any 

< ZnV -y,(f>>= 0, VV> e H N , 
or equivalently, 

TV 

£isry{x) =y]aiHi(x). (3.5) 

i=0 

3.4- Domain scaling 



It has already been mentioned in [64| that when using a spectral approach on the whole 
real line R one can possibly increase the accuracy of the computation by a suitable scaling 
of the underlying time variable t. For example, if y denotes a solution of the ordinary 
differential equation, then the rescaled function is y(t) = y{\), where I is constant. Domain 
scaling is used in several ap plications presented in next section. For more details we refer 



the interested reader to 65 



4. Applications 

In this section we apply Hermite functions collocation (HFC) method for the computation 
of Lane-Emden type equations based on the transformed Hermite functions. In general the 
Lane-Emden type equations are formulated as 

OL 

y'\x) + -y(x) + f{x)g{y) = h{x), ax > 0, (4.1) 
x 

with initial conditions 

a. ?/(0) = A, 

b. y'(0) = B, (4.2) 

10 



where a, A and B are real constants and f(x), g(y) and h(x) are some given functions. 
For other special forms of g(y), the well-known Lane-Emden equations were used to model 
several phenomena in mathematical physics and astrophysics such as the theory of stellar 
structure, the thermal behavior of a spherical cloud of gas, isothermal gas spheres and the 



theory of thermionic currents [57J, [66 



In this section we apply the Hermite functions collocation method to solve some well- 
known Lane-Emden type equations for various f(x), g{y), A and B, in two cases homoge- 
neous (h(x) = 0) and non-homogeneous (h(x) ^ 0). 

The Hermite functions are not differentiable at the point x — 0, therefore to satisfy the 
second boundary condition we can multiply the operator (13. 5p by x and add it with Bx, and 
also for satisfying the first boundary condition we add it with A as follows: 

€nV(?) = A + Bx + xi N y(x). 

Now for boundary conditions we have £,Ny(%) — A and ^£,Ny(x) = B when x tends to zero. 
To apply the collocation method, we construct the residual function by substituting y(x) by 
£,Ny(x) in the Lane-Emden Eq. (I4.ip : 

d 2 - ad- - 

Res(x) = —£ N y(x) + -— £ N y(x) + f(x)g(^ny(x)) - h(x). 

The equations for obtaining the coefficient a^s arise from equalizing Res(x) to zero at iV + 1 
transformed Hermite-Gauss points by Eq. (13.31) . i.e, transformed roots of Hermite Hn+i(x): 



Res{xj) = 0, j = 0,1,2, ...,N. (4.3) 

Solving this set of equations we have the approximating function ^y(x). We note that 
these N + 1 equations generate a set of N + 1 nonlinear equations which can be solved by 
a well-known method such as the Newton method for unknown coefficients OjS. 

4-1. The homogeneous Lane-Emden type equations 
4-1.1. Example 1 (The standard Lane-Emden equation) 

For f(x) = 1, g(y) = y m , A = 1 and B = 0, Eq. (14.11) is the standard Lane-Emden 
equation that was used to model the thermal behavior of a spherical cloud of gas acting under 
the mutual attraction of its molecules and subject to the classical laws of thermodynamics 



26, 67 



y"{x) + -y\x) + y m {x) = 0, x > 0, (4.4) 

x 

subject to the boundary conditions 

y(o) = i, 

y'(0) = 0, 
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where m > is constant. Substituting m = 0, 1 and 5 into Eq. (I4.4p leads to the exact 
solution 



1 ' f '\ / 2 \ — 

y (x) = 1 - — x , y(x) = — - — and y(x) = II + — 

respectively. In other cases there aren't any analytic exact solutions. Therefore, we apply 
the HFC method to solve the standard Lane-Emden Eq. (14.41) for m = 1.5, 2, 2.5, 3 and 4. 
To this way now we can construct the residual function as follows: 

Re Sl (x) = -f^ N y{x/l) + -^ N y(x/l) + (£ N y(x/l)) m , 

where / is a constant that is already defined in the domain scaling description. As said 
before, to obtain the coefficients a^s, Resi{x) is equalized to zero at iV + 1 transformed 
Hermite-Gauss points by Eq. (13. 3ft : 

Re Sl { Xj ) = 0, j = 0, 1, 2, N. 

By solving this set of equations, we can find the approximating function ^y(x). 

Table [1] shows the comparison of the first zeros of the standard Lane-Emden equations, 



from the present method and exact values given by Horedt [68| for m = 1.5, 2, 2.5, 3 and 4, 
respectively. 

Tables [2, [3] show the approximations of y(x) for the standard Lane-Emden for m = 3,4 
resp ectively obtained by the method proposed in this paper and those obtained by Horedt 



Table H] represents the coefficients of the Hermite functions obtained by the present 
method for m = 2, 3 and 4 of the standard Lane-Emden equation. The resulting graph of 
the standard Lane-Emden equation for m = 1.5, 2, 2.5, 3 and 4 is shown in Figured! The 
logarithmic graph of absolute coefficients of Hermite functions of standard Lane-Emden for 
m = 3 is shown in Fig. |2] and the coefficients in Table H] show that the new method has an 
appropriate convergence rate. 

4-1.2. Example 2. (The isothermal gas spheres equation) 

For f(x) = 1, g(y) = e y , A = and B = 0, Eq. (14. ip is the isothermal gas sphere 
equation. 

y"(x) + -y'(x) + e y{x) = 0, x > 0, (4.5) 

x 

subject to the boundary conditions 

y(0) = 0, 
y'(0) = 0. 

This model can be used to view the isothermal gas spheres, where the temperature remains 



constant. For a thorough discussion of the formulation of Eq. (14.51) . see 67 
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A series solution obtained by Wazwaz 29J, Liao 33J, Singh et al. 49] and Ramos 37J by 
using ADM, ADM, MHAM and series expansion respectively: 



y(x) 



1 2 a 1 4 

— x H -x 

6 5.4! 



:X° + 



122 



61.67 



-x 10 (A 6) 

21.6! ' 81.8! 495.10! ' y ' ' 

We recall that this equation has been solved by 0, 0, M, 5| with HPM, series solutions, 
ADM and HAM methods respectively. We intend to apply the HFC method to solve the 
isothermal gas spheres Eq. (14. 5 j) too. Therefore, we construct the residual function as 
follows: 

d 2 



Resi(x) 



2 d 



£ N y(x/l) + -—£ N y(x/l) + e «"»C*/0), 



dx 2 ~* 1 " i ' x 1 ' ' x dx 



where / is a constant that is already defined in domain scaling description. As said before, to 
obtain the coefficients a,s, Resi(x) is equalized to zero at iV + 1 transformed Hermite-Gauss 
points by Eq. (13. 3ft : 



Resi(xj) = 0, 



0,l,2,...,iV. 



By solving the set of equations, we have the approximating function Cnv{x). 

Tables [5] shows the comparison of y(x) obtained by the method proposed in this paper 



with (N = 30, I = 2 and k = 2) and those obtained by Wazwaz 29 



The resulting graph of the isothermal gas spheres equation in comparison to the presented 
method and those obtained by Wazwaz 29j are shown in Figure [3j The logarithmic graph 
of the absolute coefficients of Hermite functions of the standard isothermal gas spheres is 
shown in Figure HI This graph shows that the new method has an appropriate convergence 
rate. 

4.1.3. Example 3. 

For f(x) = 1, g(y) = sinh(?/), A = 1 and B = 0, Eq. (14. ip will be one of the Lane-Emden 
type equations that is we want to solve: 



y"(x) H — y\x) + sinh(y) = 0, x > 0, 
subject to the boundary conditions 

y(o) = 1, 

y'(0) = 0. 



(4.7) 



A series solution obtained by Wazwaz 29j by using Adomian Decomposition Method (ADM) 



is: 



y{x) 



(e 2 - 1) x 2 1 



l)x 4 



1 (2 e 6 + 3 e 2 - 3 e 4 - 2) x 6 



12e 



480 



30240 



+ 



1 (61 e 8 - 104 e 6 + 104 e 2 - 61) x s 



26127360 
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We intend to apply HFC method to solve equation Eq. (14.71) . 
Therefore, we construct the residual function as follows: 

d 2 - 2d- ~ 

Resi(x) = —£ N y(x/l) + ~—£ N y(x/l) + sinh (£ N y(x/l)), 
diOc iC (XiijC 

where I is a constant that is already defined in domain scaling description. To obtain the 
coefficients a^s, Resi(x) is equalized to zero at N + 1 transformed Hermite-Gauss points by 
Eq. (USD: 

Resi{xj) = 0, j = 0, 1, 2, N. 

By solving this set of equations, we have the approximating function ^y(x). 

Tables [6] shows the comparison of y(x) obtained by the new method proposed in this 
paper with (n = 10, k = 1 and 1 = 2), and those obtained by Wazwaz [29]. The resulting 
graph of Eq. (14. 7p in comparison to the presented method and those obtained by Wazwaz 



29] are shown in Figure |5J 



4- 1-4- Example 4- 

For f(x) = 1, g(y) = sin(y), A = 1 and -8 = 0, Eq. (14. ip will be one of the Lane-Emden 
type equations that is we would like to solve: 

2 

y"(x) + -y\x) + sin(y) = 0, x > 0, (4.8) 

x 

subject to the boundary conditions 

y(o) = h 

y'(0) = 0. 



A series solution obtained by Wazwaz [29J by using ADM is: 

11 11 

y(x) ~1 kix 2 H kikox 4 + fci( k 2 k1)x & 

yy ' 6 120 v 3024 1 5040 2J 

113 1 

+ hk 2 ( k 2 H k 2 )x 8 

v 3265920 1 362880 21 

v 898128000 1 2 399168000 2 23950080 1J 
where ki = sin(l) and k 2 = cos(l). 

We intend to apply the HFC method to solve this type equation (Eq. I4.8p . 
Therefore, we construct the residual function as follows: 

d 2 ~ 2 d - - 

Resi{x) = -^^nv{x/1) + ~^nv(x/1) + sin (£ N y(x/l)), 
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where / is a constant that is already defined in domain scaling description. As said before, to 
obtain the coefficients a«s, Resi(x) is equalized to zero at iV+ 1 transformed Hermite-Gauss 
points by Eq. (13. 3ft : 

Res l {x j ) = 0, j = 0, 1, 2, N. 

By solving this set of equations, we have the approximating function ^y(x). 

Tables [7J shows the comparison of y(x) obtained by the method proposed in this paper 
with (n — 15, k — 1 and 1 = 2), and those obtained by Wazwaz 29J. In order to compare 
the present method with those obtained by Wazwaz [29() the resulting graph of Eq. (14. 8 p is 
shown in Figure El 

4-1.5. Example 5. 

For f(x) = 1, g(y) = A(2e v + e y/2 ), A = and B = 0, Eq. flO) will be one of the 
Lane-Emden type equations that is to solve. 

y"(x) + -y\x) + 4(2e y + e v/2 ) = 0, x > 0, (4.9) 

x 

subject to the boundary conditions 

y(o) = o, 

y'(0) = 0. 

which has the following analytical solution: 

y(x) = -2ln(l + x 2 ). (4.10) 



This type of equation has been solved by [47|, iTOj with VIM and HPM methods respectively. 

We applied the HFC method to solve equation Eq. (14. 9p . Therefore, we construct the 
residual function as follows: 

Re Sl (x) = ^NV(x/l) + l^NV(x/l) + 4(2e^/') + e ^ x ^ 2 ), 

where I is a constant that is already defined in domain scaling description. To obtain the 
coefficients a^s, Resi(x) is equalized to zero at iV + 1 Hermite-Gauss points, i.e, roots of 
Hermite Hn+i{x): 

Res l (x j ) = 0, j = 0,l,2,...,iV. (4.11) 

By solving the set of equations, we have the approximating function ^nu(x). 

Tables [S] shows the comparison of y(x) obtained by the new method proposed in this 
paper with {n = 30, k = 2/3 and I = 2) and the analytic solution Eq. (I4.10p . In order to 
compare the present method with the analytic solution, the resulting graph of Eq. (14. 9 p is 
shown in Figure [71 
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4-1.6. Example 6. 

For f(x) = 1, g(y) = -Gy - 4yln(y), A 
Lane-Emden type equations that is: 



y"(x) + -y'(x) 
x 



Gy(x) = Ay(x) \n(y(x)), 
subject to the boundary conditions 

y(o) 



1 and B = 0, Eq. f H~Tj) will be one of the 
x > 0, (4.12) 



y'(o) 



o. 



which has the following analytical solution: 

y(x) = e x \ 



(4.13) 



This type of equation has been solved by |36l . 147] with VIM and linearization methods 
respectively. Now we apply the HFC method to solve this type of equation but in this 
model we have y(x)ln(y(x)) term that increases the order of calculation; therefore, we can 
use the transform y(x) = e z ^ in which z(x) is unknown; where upon transformed form of 
the model will become as follows: 



6 = 4z{x), x>0, 



z"(x) + (z'(x)f + -z'ix) 

X 

with the boundary conditions 

z(0) = 0, 
z'(0) = 0. 

We applied the HFC method to solve this type of equation (i. e. Eq. (I4.14p ). 
Therefore, we construct the residual function as follows: 



(4.14) 



d 2 - / d - 

Resi(x) = —£ N z(x/l) + ( —£ N z{x/l) 



+ \^nz{x/1) - 4(6v^(x/0) - 6, 



where / is a constant that is defined in domain scaling description, before. To obtain the 
coefficients a^s, Resi(x) is equalized to zero at iV + 1 transformed Hermite-Gauss points by 
Eq. (J3D: 

Res i (x j ) = 0, j = 0, 1, 2, N. 

By solving this set of equations, we have the approximating function £ N z(x) and also ^nu(x). 

Tables M shows the comparison of y(x) obtained by the method proposed in this paper 
with (n = 30, k = 6 and L = 2), and the analytic solution i. e. Eq. f)4.13p . The resulting 
graph of the Eq. (14.121) with presented method is shown in Figure [8] to make it easier to 
compare with the analytic solution . 

The logarithmic graph of absolute coefficients of Hermite functions of the standard 
isothermal gas spheres is shown in Figure El This graph shows that the new method has an 
appropriate convergence rate. 
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4-1.7. Example 7. 

For f(x) = -2(2x 2 + 3), g(y) = y, A = 1 and B = 0, Eq. gj) will be one of the 
Lane-Emden type equations that is 

y"{x) + -y'{x) - 2(2x 2 + 3)y = 0, x > 0, (4.15) 

x 

subject to the boundary conditions 

y(o) = l, 

y'(0) = 0, 

which has the following analytical solution: 

y{x) = e x \ (4.16) 



This type of equation has been solved by |36|, |47j, 170J with linearization, VIM and HPM 
methods respectively. We applied the HFC method to solve the equation ( I4.15p . 
Therefore, we construct the residual function as follows: 

Resi(x) = ^£ N y{x/l) + -^ N y(x/l) - 2{2x 2 + 3)£ N y(x/l), 
ax^ x ax 

where / is a constant that is already defined in domain scaling description. As said before, to 
obtain the coefficients a^s, Resi(x) is equalized to zero at iV + 1 transformed Hermite-Gauss 
points by Eq. (13. 3ft : 

Res l (x j ) = 0, j = 0, 1, 2, N. 

By solving this set of equations, we have the approximating function ^Ny(x). 

Tables [TO] shows the comparison of y(x) obtained by the new method proposed in this 
paper with (n = 30, k = 6 and 1 = 2), and analytic solution Eq. (I4.16p . 

The resulting graph of Eq. (14.151) is shown in Figure [101 Thus it will be easy for the 
reader to compare the new result with the exact solution. The logarithmic graph of absolute 
coefficients of Hermite functions of Eq. (I4.15P is shown in Figure [TTJ This graph shows that 
the current method has an appropriate convergence rate. 

4-2. The non-homogeneous Lane-Emden type equations 
4-2.1. Example 8. 

For f{x) = x, g{y) = y, h(x) = x 5 - x 4 + 44x 2 - 30x, a = 8, A = and B = 0, Eq. <KT\f 
will be one of the Lane-Emden type equations that is absorbing to solve, 
g 

y"(x) + -y'(x) + xyix) = x 5 - x 4 + 44x 2 - 30x, x > 0, (4.17) 

x 

subject to the boundary conditions 

y(0) = 0, 
y'(0) = 0, 
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which has the following analytical solution: 

y (x) = x 4 - a; 3 . (4.18) 

This type of equation has been solved by 36|, |42|, |7l|, lZ2| with linearization, HPM, HAM and 
TSADM methods respectively. We applied the HFC method to solve this type equation Eq. 
( I4.17p . Therefore, we construct the residual function as follows: 

d 2 — 8 d — — 

Resi(x) = —£ N y(x/l) + -—£ N y(x/l) + x£ N y{x/l) - x 5 + x 4 - AAx 2 + 30a;, 
dx A x ax 

where I is a constant that is already defined in domain scaling description. To obtain the 
coefficients a^s, Resi(x) is equalized to zero at N + 1 transformed Hermite-Gauss points by 
Eq. (J33): 

Res l (x j ) = 0, j = 0, 1, 2, N. 

By solving this set of equations, we have the approximating function £ N y(x). 

Tables [TT] shows the comparison of y(x) obtained by the method proposed in this paper 
with (n = 30, k = 2/3 and 1 — 2), and analytic solution Eq. (14.18p . The resulting graph of 
the Eq. (14.1 7p for the presented method and the analytic solution are shown in Figure [12j 

4-2.2. Example 9. 

For f(x) = 1, g(y) = y, h(x) = 6 + 12x + x 2 + x 3 , A = and B = 0, Eq. P~U) will be 
one of the Lane-Emden type equations that is 
2 

y"(x) + -y'ix) + y(x) = 6 + 12x + x 2 + x 3 , x>0, (4.19) 
x 

subject to the boundary conditions 

2/(0) = 0, 
y'(0) = 0, 

which has the following analytical solution: 

y(x)=x 2 + x 3 . (4.20) 

This equation has been solved by 36|, |47|, |70|, |72| with with linearization, VIM, HPM and 
TSADM methods respectively. We applied the HFC method to solve Eq. f)4.19p . 
Therefore, we construct the residual function as follows: 

d 2 - 2d- ~ 

Resi(x) = —£, N y(x/l) + -—£ N y(x/l) + £ N y{x/l) - 6 - 12x - x 2 - x 3 , 
dx z x dx 

where / is a constant that is already defined in domain scaling description. As said before, to 
obtain the coefficients a«s, Resi(x) is equalized to zero at iV+ 1 transformed Hermite-Gauss 
points by Eq. (13. 3ft : 

Res l {x j ) = 0, j = 0, 1, 2, N. 

By solving this set of equations, we have the approximating function ^y(x). Tables [T21 
shows the comparison of y(x) obtained by the method proposed in this paper with (n = 30, 
k = 2/3 and / = 2), and analytic solution Eq. (I4.20p . The resulting graph of Eq. (I4.19P 
with presented method in comparison to the analytic solution is shown in Figure [13j 
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5. Conclusions 



The Lane-Emden equations describe a variety of phenomena in theoretical physics and 
astrophysics, including the aspects of stellar structure, the thermal history of a spherical 



cloud of gas, isothermal gas spheres, and thermionic currents [57 



have been considered by many mathematicians as mentioned before [55 



Lane-Emden equations 
The fundamental 



goal of this paper has been to construct an approximation to the solution of nonlinear Lane- 
Emden type equations in a semi-infinite interval. A set of Hermite functions is proposed to 
provide an effective but simple way to improve the convergence of the solution by the collo- 
cation method. The validity of the method is based on the assumption that it converges by 
increasing the number of collocation points. A comparison is made among the exact solution 



68 and the series solutions of Wazwaz 



Liao 33 



and the numerical solutions of Horedt 

Singh et al. 49j and Ramos 37j and the current work. It has been shown that the present 
work provides acceptable approach for Lane-Emden type equations. Also it was confirmed 
by logarithmic figures of absolute coefficients, this approach has exponentially convergence 
rate. In total an important concern of spectral methods is the choice of basis functions; 
the basis functions have three different properties: easy to compute, rapid convergence and 
completeness, which means that any solution can be represented to arbitrarily high accuracy 
by taking the truncation N to be sufficiently large. 
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Table 1: Comparison of the first zeros of standard Lane-Emden equations, for the present method and exact 
numerical values given by Horedt [681 ] 

m N k I Present method Exact value 



1.5 


4 


1 


3.74224350 


3.65375374 


3.65375374 


2 


10 


1 


1.97027600 


4.35287460 


4.35287460 


2.5 


10 


1 


1.97668316 


5.35527546 


5.35527546 


3 


20 


1 


1.86927585 


6.89684862 


6.89684862 


4 


12 


1/3 


1.97137830 


14.9715463 


14.9715463 



Table 2 


: Comparison of y(x) values of standard Lane-Emden equation, 


for the present method and exact 


values g 


,iven by Horedt I68j. for m=3 






X 


Present method 


Exact value 


Error 


0.0 


1.00000000 


1.0000000 


O.OOe + 00 


0.1 


0.99833720 


0.9983358 


1.40e - 06 


0.5 


0.95984209 


0.9598391 


2.99e - 06 


1.0 


0.85505959 


0.8550576 


1.99e-06 


5.0 


0.11082019 


0.1108198 


3.89e - 07 


6.0 


0.04373912 


0.0437380 


1.12e - 06 


6.8 


0.00417826 


0.0041678 


1.05e-05 


6.896 


0.00003610 


0.0000360 


9.79e - 08 


Table 3 


: Comparison of y(x) values of standard Lane-Emden equation, 


for the present method and exact 


values g 


iven by Horedt |68j, for m=4 






X 


Present method 


Exact value 


Error 


0.0 


1.0000000 


1.0000000 


O.OOe + 00 


0.1 


0.9985876 


0.9983367 


2.51e-04 


0.2 


0.9936339 


0.9933862 


2.48e - 04 


0.5 


0.9605160 


0.9603109 


2.05e - 04 


1.0 


0.8610072 


0.8608138 


1.93e-04 


5.0 


0.2358368 


0.2359227 


8.59e-05 


10.0 


0.0596105 


0.0596727 


6.22e-05 


14.0 


0.0083058 


0.0083305 


2.47e - 05 


14.9 


0.0005759 


0.0005764 


4.59e - 07 
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Table 4: Coefficients of the Hermite functions of the standard Lane-Emden equations for m — 2,3 and 4 



respectively 



i a j 




i a i 




m = 2 




rn = 3 




m — 4 








m = 2 


m = 3 




m — 4 





-5.2841135322e 


- 01 


-4. 4099373672s 


- 01 


-3.8511246127e 


- 01 




13 




-2.4548768173e 


- 02 




1 


-2.0672313847e 


- 01 


-1.5728415017e 


- 01 


1.16850585566 


- 01 




14 




-1.79164202816 


- 02 




2 


-2.1013493211e 


- 01 


-1.7607131187e 


- 01 


-1.65766227136 


- 01 




15 




-1.02002272586 


- 02 




3 


-1.2898718939e 


- 01 


-1.1378421470e 


- 01 


9.2854106306e 


- 03 




16 




-6.5268030714e 


- 03 




4 


-1.3634530855e 


- 01 


-1.2995159559e 


- 01 


-7.1551541010e 


- 02 




17 




-2.7962896018e 


- 03 




5 


-8.7619773995e 


- 02 


-9. 6296863459s 


- 02 


-9.8809929827e 


- 03 




18 




-1.5765572392e 


- 03 




6 


-7.2750465809e 


- 02 


-9.8373526479e 


- 02 


-4.8372346356e 


- 02 




19 




-3.7895054857e 


- 04 




7 


-3.9156883681e 


- 02 


-7.9430021072e 


- 02 


-9.4556001733e 


- 03 




20 




-2.4542997154e 


- 04 




8 


-2.6813942695e 


- 02 


-7.8340439572e 


- 02 


-2.6810185942e 


- 02 














9 


-9.5249929620e 


- 03 


-6. 2915940155a 


- 02 


-6.6016540826e 


- 03 














10 


-4.1991804282e 


- 03 


-5.7157720774e 


- 02 


-1.19100532776 


- 02 














11 






-4.3579589433e 


- 02 


-1.1402951223e 


- 03 














12 






-3.6177724390e 


- 02 


-3.1134305650e 


- 03 















Table 5: Comparison of y(x) , between present method and series solution 


^iven by Wazwaz [29( for isother- 


mal gas sphere equation 






x Present method 


Wazwaz 


Error 


0.0 0.0000000000 


0.0000000000 


0.00e + 00 


0.1 -0.0016664188 


-0.0016658339 


5.85e - 07 


0.2 -0.0066539713 


-0.0066533671 


6.04e - 07 


0.5 -0.0411545150 


-0.0411539568 


5.58e-07 


1.0 -0.1588281737 


-0.1588273537 


8.20e-07 


1.5 -0.3380198308 


-0.3380131103 


6.72e - 06 


2.0 -0.5598233120 


-0.5599626601 


1.39e-04 


2.5 -0.8063410846 


-0.8100196713 


3.68e-03 



Table 6: 
3 


Comparison of y{x) , between present method and series solution 


jiven by Wazwaz [29| for Example. 


X 


Present method 


Wazwaz 


Error 


0.0 


1.0000000000 


1.0000000000 


0.00e + 00 


0.1 


0.9981138095 


0.9980428414 


7.10e-05 


0.2 


0.9922758837 


0.9921894348 


8.64e - 05 


0.5 


0.9520376245 


0.9519611019 


7.65e - 05 


1.0 


0.8183047481 


0.8182516669 


5.31e-05 


1.5 


0.6254886192 


0.6258916077 


4.03e - 04 


2.0 


0.4066479695 


0.4136691039 


7.02e - 03 
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Table 7: Comparison of y(x) , between present method and series solution given by Wazwaz [29( for Example. 



4 



X 


Present method 


Wazwaz 


Error 


0.0 


1.0000000000 


1.0000000000 


0.00e + 00 


0.1 


9986051425 


0.9985979358 


7.21e - 06 


0.2 


n 9944062706 

\J • U J1 lU \J £-1 \ \J\J 


0.9943962733 


l.OOe - 05 


0.5 


9651881683 


0.9651777886 


1.04e - 05 


1.0 


8636881301 


0.8636811027 


7.03e - 06 


1 5 

-L . "J 


70505241 03 


0.7050419247 


1.05e - 05 


2 


5064687568 


0.5063720330 


9.67e - 05 




Table 8: Comparison of y(x) , 


between present method and exact solution for Example. 5 


X 


Present method 


Exact value 


Error 


0.00 


0.0000000000 


0.0000000000 


O.OOe + 00 


0.01 


-0.0001970587 


-0.0001999900 


2.93e - 06 


0.10 


-0.0198967225 


-0.0199006617 


3.94e - 06 


0.50 


-0.4462840851 


-0.4462871026 


3.02e - 06 


1.00 


-1.3862934297 


-1.3862943611 


9.31e-07 


2.00 


-3.2188763248 


-3.2188758249 


5.00e-07 


3.00 


-4 6051709964 


-4.6051701860 


8.10e-07 


4.00 


-5.6664274573 


-5.6664266881 


7.69e - 07 


5.00 


-6 5161937402 


-6.5161930760 


6.64e - 07 


6.00 


—7 2218363729 


-7.2218358253 


5.48e - 07 


7.00 


—7 8240461812 


-7.8240460109 


1.70e - 07 


8.00 


—8 3487734467 


-8.3487745398 


1.09e - 06 


9 00 


—8 8134506165 


-8.8134384945 


1.21e-05 


10 00 

X \J • \J\J 


—9 2302027821 


-9.2302410337 


3.83e-05 




Table 9: Comparison of y(a;) , 


between present method and exact solution for Example. 6 


X 


Present method 


Exact value 


Error 


0.00 


1.0000000000 


1.0000000000 


O.OOe + 00 


0.01 


1.0000999826 


1.0001000050 


2.24e - 08 


0.02 


1.0004000642 


1.0004000800 


1.58e-08 


0.05 


1.0025031064 


1.0025031276 


2.12e-08 


0.10 


1.0100501492 


1.0100501671 


1.79e-08 


0.20 


1.0408107527 


1.0408107742 


2.15e-08 


0.50 


1.2840253862 


1.2840254167 


3.05e-08 


0.70 


1.6323161777 


1.6323162200 


4.23e - 08 


0.80 


1.8964808279 


1.8964808793 


5.14e-08 


0.90 


2.2479078937 


2.2479079867 


9.29e - 08 


1.00 


2.7182819166 


2.7182818285 


8.81e-08 
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Tabic lOi Comparison of 


, between present method and exact solution for Example. 7 


X 


Present method 


Exact value 


Error 


0.00 


1.0000000000 


1.0000000000 


O.OOe + 00 


0.01 


1.0000999826 


1.0001000050 


2.24e - 08 


0.02 


1.0004000642 


1.0004000800 


1.58e-08 


0.05 


1.0025031065 


1.0025031276 


2.12e-08 


0.10 


1.0100501493 


1.0100501671 


1.78e-08 


0.20 


1.0408107533 


1.0408107742 


2.09e-08 


0.50 


1.2840253904 


1.2840254167 


2.62e-08 


0.70 


1.6323161872 


1.6323162200 


3.27e - 08 


0.80 


1.8964808414 


1.8964808793 


3.79e - 08 


0.90 


2.2479079319 


2.2479079867 


5.48e - 08 


1.00 


2.7182818260 


2.7182818285 


2.51e-09 



Table 11: Comparison of y(x) , between present method and exact solution for Example. 8 


X 


Present method 


Exact value 


Error 


0.00 


0.0000000000 


0.0000000000 


O.OOe + 00 


0.01 


-0.0000009321 


-0.0000009900 


5.79e - 08 


0.10 


-0.0009008409 


-0.0009000000 


8.41e-07 


0.50 


-0.0625021958 


-0.0625000000 


2.20e - 06 


1.00 


-0.0000008284 


0.0000000000 


8.28e-07 


2.00 


8.0000001732 


8.0000000000 


1.73e-07 


3.00 


54.0000002074 


54.0000000000 


2.07e - 07 


4.00 


192.0000000368 


192.0000000000 


3.68e - 08 


5.00 


499.9999998091 


500.0000000000 


1.91e-07 


6.00 


1079.9999995264 


1080.0000000000 


4.74e - 07 


7.00 


2058.0000004141 


2058.0000000000 


4.14e-07 


8.00 


3584.0000093640 


3584.0000000000 


9.36e-06 


9.00 


5831.9999560359 


5832.0000000000 


4.40e - 05 


10.00 


8999.9996608001 


9000.0000000000 


3.39e - 04 
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Table 12: Comparison of y{x) , between present method and exact solution for Example. 9 



X 


Present method 


Exact value 


Error 


0.00 


0.0000000000 


0.0000000000 


O.OOe + 00 


0.01 


0.0000995275 


0.0001010000 


1.47e - 06 


0.10 


0.0109981790 


0.0110000000 


1.82e-06 


0.50 


0.3749985918 


0.3750000000 


1.41e-06 


1.00 


1.9999987524 


2.0000000000 


1.25e-06 


2.00 


11.9999993068 


12.0000000000 


6.93e-07 


3.00 


35.9999999242 


36.0000000000 


7.58e - 08 


4.00 


80.0000003071 


80.0000000000 


3.07e - 07 


5.00 


150.0000003207 


150.0000000000 


3.21e-07 


6.00 


252.0000000974 


252.0000000000 


9.74e - 08 


7.00 


391.9999997951 


392.0000000000 


2.05e-07 


8.00 


575.9999992644 


576.0000000000 


7.36e - 07 


9.00 


810.0000046092 


810.0000000000 


4.61e-06 


10.00 


1099.9999875537 


1100.0000000000 


1.24e - 05 
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Figure 1: Graph of standard Lane-Emden equation for m = 1.5, 2, 2.5, 3 and 4 
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Figure 4: Logarithmic graph of absolute coefficients \at \ of Hermite function of isothermal gas sphere equation 
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[ HFC ♦ Wazwaz | 

Figure 5: Graph of equation example 3 in comparing the presented method and Wazwaz solution [2' 
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| HFC ■ Exact | 

Figure 7: Graph of equation example 5 in comparing the presented method and analytic solution 




X 

I HFC ■ Exact] 

Figure 8: Graph of equation example 6 in comparing the presented method and analytic solution 
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Figure 9: Logarithmic graph of absolute coefficients |dj| of Hermite function of example 6 
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Figure f 0: Graph of equation example 7 in comparing the presented method and analytic solution 
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Figure 11: Logarithmic graph of absolute coefficients \<n\ of Hermite function of equation of example 7 
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Figure 12: Graph of equation example 8 in comparing the presented method and analytic solution 
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HFC ■ Exact 



Figure 13: Graph of equation example 9 in comparing the presented method and analytic solutio 
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